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By Bernard Budiansky and Edwin T. Kruszewski 
SUMMARY 


The variational principle, differential equations, and boundary con- 
ditions considered appropriate to the 'analysis of transverse vibrations 
of hollow thin-walled cylindrical beams are shown. General solutions 
for the modes and frequencies of cantilever and free-free cylindrical 
beams of arbitrary cross section but of uniform thickness are given. 

The combined influence of the secondary effects of transverse shear 
deformation, shear lag, and longitudinal inertia is shown in the form 
of curves for cylinders of rectangular cross section and uniform thick- 
ness. The contribution of each of the secondary effects to the total 
reduction in the actual frequency is also indicated. 


INTRODUCTION 


The elementary theory of bending vibration is often inadequate for 
the accurate calculation of natural modes and frequencies of hollow, thin- 
walled cylindrical beams. Such secondary effects as transverse shear 
deformation, shear lag, and longitudinal inertia, which are not considered 
in the elementary theory of lateral oscillations, can have appreciable 
influence, particularly on the higher modes and frequencies of vibration. 
The effects of transverse shear deformation and of rotary (rather than 
longitudinal) inertia have been studied by many on the basis of the 
original investigations of Rayleigh (reference l) and Timoshenko (refer- 
ence 2). Anderson and Houbolt (reference 3) have presented a procedure 
for including the effects of shear lag in the numerical calculation of 
modes and frequencies of box beams of rectangular cross section. How- 
ever, there does not appear to exist a general solution for the vibra- 
tion of hollow beams that incorporates the influence of all the secondary 
effects mentioned. 

The purpose of the present paper is threefold: First, to exhibit 

the variational principle, differential equations, and boundary condi- 
tions appropriate for the analysis of the uncoupled bending vibration 
of hollow thin-walled cylindrical beams; second, to give general solutions 
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for cantilever and free-free cylinders of arbitrary cross section but of 
uniform thickness; and finally, to show quantitatively the influence of 
the secondary effects by means of numerical results for hollow beams of 
rectangular cross section of various lengths, widths, and depths. 
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SYMBOIS 


cross-sectional area 
Fourier coefficient 

effective shear-carrying area ’ - 

parameter defined in equation (30) 

constant 

modulus of elasticity 
shear modulus of elasticity 
moment of inertia 

geometrical parameter defined in equation (29) 

length of cantilever beam, half-length of free-free beam 

parameter defined in equation (38) 

maximum kinetic energy 

maximum strain energy 

half-depth of rectangular beam 

half-width of rectangular beam 

Fourier series coefficients 


integers 


frequency coefficient 
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coefficient of shear rigidity 


coefficient of rotary inertia 



perimeter of cross section 

distance along periphery of cross section (see fig. l) 
wall thickness 

longitudinal displacement in x-direction 
vertical displacement in y-direction 
longitudinal coordinate 
vertical coordinate 

y-coordinate of center of gravity of cross section 
shear strain 


longitudinal strain 

inclination of normal with vertical (see fig. l) 

Lagrangian multiplier 

mass of beam per unit length 

mass density of beam 

longitudinal direct stress 

shear stress 

natural frequency of beam 

natural frequency of beam calculated from elementary beam 
theory 



it 


NACA TN 2682 


BASIC EQUATIONS 


Assumptions .- The problem to be considered is that of the natural 
bending vibration of a thin-walled hollow cylindrical beam whose cross 
section is symmetrical, about at least one axis (see fig. l). The 
transverse vibration is supposed to take place in the direction of this 
axis of symmetry of the cross section so that no torsional oscillations 
are induced. 

In the present analysis, the following simplifications are introduced: 

(a) Changes in the size and shape of the cross section are neglected 

(b) Stress and strain are assumed to be uniform across the wall 
thickness 

(c) The small effect of circumferential stress upon longitudinal 
strain is neglected 

In accordance with statements (a) and (b), the distortions of the 
vibrating beam are completely described by the vertical displacement w(x) 
of a cross section and the longitudinal displacement u(x,s) of each 
point of the median line of the beam wall. 

The longitudinal and shear strains are given in terms of u(x,s) 
and w(x) as 


and 


e = 

Cx 8x 


7 = + dw sin e 

xs 8s dx 


( 1 ) 


( 2 ) 


and the corresponding stresses become 


<J X = E — 
• X dx 

and 


(3) 


T *s = f sm ») W 

where 9 is the inclination of the normal with the vertical (see fig. l), 
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In elementary beam theory, where the effects of all shear distortion 
are neglected, the longitudinal distortion u(x,s) is related to the 
vertical displacement w(x) by - 

u(x,s) = (y - y)—£ 
dx 

where y is the y-coordinate of the center of gravity of the' cross sec- 
tion. In the present paper, however, u(x,s) is allowed to be perfectly 
general, so that shear distortions (and consequently the so-called shear- 
lag and transverse-shear-deformation effects) are fully taken into account. 
Furthermore, because cross sections are not constrained to remain plane, 
the inertia effect associated with motion in the longitudinal direction is 
more properly designated as the effect of longitudinal inertia rather than 
rotary inertia. 

Variational principle and geometrical boundary conditions .- The varia- 
tional equation to be written is appropriate to beams whose ends are either 
fixed, simply supported, or free. For some such beam vibrating in a 
natural mode, the maximum strain energy is 

u "t X/K^) 2 tdstot ii o ^ G ® t S slne ) 2t<i0d " (5) 

where u(x,s) and w(x) are the amplitudes of displacement for the par- 
ticular mode considered. The maximum kinetic energy is 


T = i J (jf pt<u 2 w 2 ds dx + i J (jl ptCD 2 u 2 ds dx (6) 

where cjd is the natural frequency of the mode under consideration and 
p is the mass density of the beam. The second term in equation (6) con- 
stitutes the contribution of longitudinal inertia to the kinetic energy. 

Then a natural mode of vibration must satisfy the variational 
equation 

5(U - T) = 0 (7) 


where the variation is taken independently with respect to u(x,s) and 
w(x) and with the provision that both u(x,s) and w(x) 'must satisfy 



the geometrical boundary conditions of the problem; furthermore, u(x,s) 
must be periodic in s with a period equal to the perimeter p. The 
geometrical boundary conditions are w = 0 and u = 0 at a fixed end 
and just w = 0 at a simply supported end. At a free end no geometrical 
boundary conditions are imposed. 


Differential equations and natural boundary conditions .- Equations (5) 
( 6 ), and. ( 7 ) in conjunction with the -usual procedure of the calculus. of 
variations yield the following simultaneous integrodifferential equations 
for u and w: 


Et 




+ — sin 6 
dx 



= Q 


( 8 ) 


^ Gt^— ^ sin 6 + sin^ejds + po)^w = 0 


( 9 ) 


where 


p = (i? pt ds 


( 10 ) 


and the boundary equations at each end of the beam are 



Su ds 


0 


( 11 ) 




dw 

dx 


sin 6 1 sin 6 ds 6w = 0 


( 12 ) 


At a fixed end, both boundary equations (ll) and (12) are satisfied 
by virtue of the fact that the geometrical boundary conditions require 
that both 6u and Sw be zero. At a simply supported end Sw = 0, but, 
since 6u(x;s) is perfectly arbitrary, the variational process forces 
the equality 
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Finally, at a free end, since there are no geometrical constraints, both 
Su and 5w are arbitrary and hence the variational process forces, in 
addition to equation (13), the equality 

^ + — sin 0^sin 6 ds = 0 (1*0 

Equations (13) and (l4) constitute so-called "natural boundary 
conditions" because they are automatically satisfied as the result of 
a variational process. Equation (13) is recognized as the condition of 
zero longitudinal direct stress while equation (l4) simply stipulates 
that the total vertical shear force vanish. 

Thus to summarize, the appropriate boundary conditions required for 
the solution of equations (8) and (9) are 

Fixed end: 

w = 0 
u = 0 


Simply supported end: 


w = 0 


Et 0 

OX 


Free end: 



dw . \ 

— sm 01 
dx / 


sin 9 ds = 0 



0 


The integrodifferential equations (8) and (9), which specify 
equilibrium in the longitudinal and transverse directions, respectively, 
can, of course, be written directly without recourse to the variational 
principle. 
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GENERAL SOLUTIONS FOR CYLINDERS OF .UNIFORM WALL THICKNESS 


The following exact solutions for cylinders of uniform wall thick- 
ness are carried out by means of Fourier series in conjunction with the 
application of the variational condition (equation ( 7 )). This procedure, 
which does not require explicit consideration of the natural boundary 
conditions, was believed to be more expedient than a direct attack upon 
the simultaneous integrodifferential equations (8) and (9) and all their 
associated boundary conditions. 

Cantilever beam .- The geometrical boundary conditions, for a canti- 
lever beam, as previously shown, are w(0) = u(0,s) = 0 (see fig. l). 
Appropriate assumptions for the displacements w(x) and u(x,s) are 

w(x) = c + b n cos (15) 

n=l,3,5 2L 

and 

00 ( 00 

/ s V- V ' . mux 2nns f ^ \ 

u(x,s) = 2_ 2_ a mn sin cos U6) 

m=l, 3,5 n=0,l, 2 P 


The condition u( 0 ,s) = 0 is satisfied by each term of equation (l6) 
the condition 


w( 0 ) = C + b n = 

n=l,3,5 


0 


(17) 


is introduced into the variational procedure by means of the Lagrangian 
multiplier method. The choice of the particular trigonometric functions 
used in the Fourier series ( 15 ) and (l6) was guided by consideration of 
the orthogonality required for the simplification of expressions in the 
strain energy. The constant C is needed in the expression for w(x) 
in order that w(L) be unrestricted. 
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Using equations (15) and (l6) in equations (5) and (6) yields 


u - T -i r$A £- £ a »I cos lr 

2J 0 J \m=l,3,5 n=0,l,2 2L 2L 


2nns \ 2 

cos ds dx + 


1 / V «,( £ £ — a* 

2J 0 J \tt=l,3,5 n=0,l,2 p 


. mux . 2nns , 

sm sm + 

2L p 


sin 9 ) -b_ sin S«x f ds dx - 

- n 2L 2L ' 


n=l,3,5 


I 


H.L *• ” 


ds dx 


Y f $ M Ti II ^ sin It 1 cos ^r)* 

2 J 0 J \m=l,3,5 n=0 ,1,2 L P / 


ds dx (18) 


To make equation (l8) stationary and at the same time satisfy the con- 
straining relationship 


<P 


♦ r 

n=l73,5 . 


b n = 


(19) 


it is sufficient to set 


5(U - T - Xcp) .= 0 


( 20 ) 
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where the variation is with respect to the a’s, b's, and C considered 
as independent variables; here X is a Lagrangian multiplier. This varia- 
tional process results in the following equations: 


8(U - T) 
8b ^ 


1-1 

| bi - (-1) 2 || puu 2 C - X. 

= 0 

(i =a, 3, 5, . . .) (21) 



8(U - T) 
aa ij 


Et (gf T a iJ 


+ Gt 






+ Gt 


ij* £ 


A j b i 


2 , Lp 
pt -£■ 
4 


ij 


= 0 


(i = 1, 3, 5, . . 

J = 0, 1 , 2, . . .') (22) 


8(U - T) 

8c 


, 2 
- X = -\XiT 

8c 


£ - 
rr nn 


n-1 

2 


n=l,3,5 


(-1) 2 b n - pu^LC - X 


= 0 


(23) 


where 




. a . 2nits , 
sin 6 sm ds 




t sin^Q ds 


(24) 


(25) 
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With the use of the nondimensional parameters 


k B 2 = Hi! ^ 
El 


( 26 ) 


*s 2 --“-5 
A S GL 


(27) 


k 2 . _I 
■ RI 


ptL 2 AL 2 


(■28) 


and 


K‘ 


2 _ 161 


A S P 


(29) 


B i 2 - 12 - k RlV(f) 2 

equations (21), (22), and (23) may be reduced to 


(30) 


00 

V - inn 2 Lt l/i*\2 1 v 1 , 2, 

~ l,t -„2 A-. An in 2V2 J b i " ? k B^i 


n=67T,2 4ks 2 ^ “ 2 \ 2 / ^2 


i-1 

(-1) 2 2 k 2 c . A 

i« . El. 


= 0 


(i - 1, 3, 5, . . ’.) (31) 


( k S 2 Bi 


2 * 


jijbi = 0 


(i = 1, 3, 5, 

j = 0, 1, 2, 


•) (32) 


n-1 


K B 


21 J? (-!) 2 b n + k B 2 C + = 


n=l,3,5 


(33) 
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For j = 0, equation (32) becomes 


ks 


2 



' k B 2k RI 



a i0 = 0 


(i = l, 3, 5, . . .) (3b) 


Equation (3b) is not coupled to any of equations (3l) to (33). A given 
value of a iQ corresponds to the amplitude of the ith mode of longi- 
tudinal oscillation, and if this value of a^Q is not equal to 0, 

then equation (3*0 simply gives the frequency of this longitudinal mode. 
Consequently those equations in equation (32) for values of j' = 0 are 
not associated with transverse bending and so are ignored henceforth. 
For the remaining values of j (that is, j ^ 0) equation (32) yields 


1 J 


XL A J 1J 

+ K 2 j 2 


(i = l, 3, 5, . . .; 
j = 2, 3 > • . .) (35) 


Substituting the expression for a^ in equation (35) into equation (31) 
and solving for b^ gives 


where 


(-D 


i-1 
2 2 v _2n 


b-j = 


lit 


k B C- + 


L3x. 

El 


N-: 


(i = 1, 3, 5, . . .) (36) 


Ni = 


i 2 it 2 1 


£ 


i k 2 


k s 2 n=l,2, 3 ks 2 ®! 2 + K^ 2 2 B 


(37) 


In the appendix this expression for 1^ is shown to be equivalent to 


i 2 it k p i 2 « 2 A k o V A n 2 n p 

1 = — 1 ‘ IT j; ^ n .fo 3 7 kV) - 2 k B (38) 


32 
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Since the series in equation ( 38 ) is considerably more quickly convergent 
than that in equation (37)> equation ( 38 ) should be used in actual 
numerical calculations of N^. 

Substitution of equation ( 36 ) into equation (33) and the constraining 
relationship equation ( 19 ) gives the following two homogeneous equations 
in C and X: 




1 + 


" b2 rX,5 ^ \ 


c + 


1 + k 


B 


00 

z 


n=l,3,5 


n-1 

(-1) 2 — — 

n« N n 


LfX 

El 


= 0 


(39a) 


1 + k B 2 



n=l,3,5 


(-D 


n-1 

2 


2 _ 1 _ 
nn N n 


C 



(39b) 


Finally the condition for a nontrivial solution for C and X gives 
the frequency equation 


*b £ 


1 + 


OO 

k B 2 y~ 

n=l,3,5 




1 + k B 2 


OO 

Z 


n-1 

(-D 2 


2_ _1 
nn N n 


w z 

n=l,3,5 


n-1 

(-D 2 


2 _ _ 1 _ 
nn N n 


= 0 (ko) 



n=l,3,5 


1_ 

N n 


which the frequency parameter kg must satisfy. Since the terms of the 
infinite series which appear in the frequency equation contain kg 
itself, the roots of equation (4o) are most conveniently found by trial. 
Fortunately the infinite series in equation (40) as well as the series 
in the definition of converge rapidly so that only a few terms are 

needed to evaluate them with sufficient accuracy. 
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Once k B has been determined for a particular mode, the corre- 
sponding mode shape can be found by letting C = 1 and solving either 
of equations (39) for X and then finally evaluating and a ± i 

successively from equations ( 36 ) and (35). 

Free -free beam - symmetrical modes .- If the origin of a free-free 
beam of length 2L is taken at the midspan (see fig. l), the form of 
the Fourier series assumed for w(x) and u(x,s) when the beam is 
undergoing a symmetrical mode of vibration may be exactly the. same as 
that assumed for the cantilever beam of length ,L .( see equations (15) 
and (16)). The only difference in. the ensuing calculations is that the 
constraining condition (19) is not introduced. Consequently, it can be 
readily seen that the frequency equation for the symmetrically vibrating 
free-free beam is obtained from equation (39a) by setting X = 0 and 
is 




1 + 


n=^J,5 (nit) »„ 


= 0 


(41) 


After a particular root k B is found from equation (4l), the shape 

of the corresponding symmetrical free-free mode may be obtained from 
equations ( 36 ) (with X = 0) and equation (35). 

Free-free beam - antisymmetrical modes .- Consider a free -free beam 
of length 2L undergoing antisymmetrical vibrations. Explicit con- 
sideration need be given only to the right half of : the beam (see fig.- l), 
and for this half -beam the only geometrical boundary condition that must 
be imposed is that w(0) = 0. The spanwise displacement u(0,s) is 
unrestrained by virtue of antisymmetry. 

Appropriate assumptions for the displacements w(x) and u(x,s) 
are then . ■ 


w(x) 



sin 


nitx 

2L 


+ Cx 


(42) 


and • 


00 00 

■ z z 

Ok n = 1 9 


o ^ mnx 
a mn cos 2L COs 


2nTts 

P 


( 43 ) 


u(x,s) 
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The linear portion Cx of the expression for w(x) is needed in 
order to give the heam sufficient freedom at the tip (x = L). The choice 
of the particular trigonometric function in the series expansion for 
u(x,s) was, as in the case of the cantilever beam, guided by considera- 
tion of the orthogonality required for the simplification of the expres- 
sions in the strain energy. The zeroeth term in the series for u(x,s) 
in the s-direction was omitted because it only leads to the frequency 
equation for longitudinal oscillations. 

Using equations (42) and (43) in equations (5) and (6) yields 



1 

2 



Gt 


£ 


m=0,2,4 



n=l,2. 


3 


00 

= T 9 


m * 

' a mn 2L 


nutx 

sin cos 

2L 


, as to + 

P J 



mjtx 

cos sm 

2L 


2mts 

P 


+ ■ 



The variation of equation (44) with respect to the a's, b's, and 
C gives, after suitable simplification, 

(Bi 2 k s 2 + K 2 j 2 ) aiJ - K 2 £ Ajijbi = 0 ( i = 2, 4, 6, . . . ; 

J “ 2 j 3 j * • • ) ( 45 ) 
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’^ 2 - *bV :V(ff 


8 oj ' "" S A J JC ' 0 


(0 = 1 , 2 , 3 , . . .) (^) 


Y“ 1 Lt ‘inn 2 11 /i rt \2, 12, 

n ln ' W\T) bi + 2 k B b i - 


( - l)l/2 j| k B 2 CL = 0 


(i.= 2, 4, 6, . (4?) 


E :=5 


n=l,2,3 


1 Lt . CL 

ks 2 A S nnA n a On ~ 


- (-D n/ \ + 

B o~T ^ a* n 

n=2,4,6 


i k R 2 CL = 0 

3 B 

From equation (45) 


K2 fe 


1J 


B^ks 2 + K^j 


b, 

2,2 1 




( L — 2 , 4 , 6, , . • J 

j = l, 2, 3 , . . .) (49) 


which, except for sign, is the same expression as that obtained for the 
cantilever and symmetrically vibrating free-free beams (equation (35)). 
From equation (46) ‘ 


K‘ 


2 P_ 
2n 


A- 




a. 


Oj 


B 0 2 kg 2 + K 2 j 2 


(j — 'I 7 2 , 3> • • •) ( 50 ) 
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Substitution of equation (49) into equation (4-7) gives 

bi = -(-l) 1 / 2 ^ CL (i = 2, 4, 6, . . .) ( 51 ) 

where is defined in equation (.37)* 

Substitution of equations (50) and (51) into equation (48) and 
simplification gives as the frequency equation for the antisymmetrically 
vibrating free -free beam 



( 2\2 2 _ A. 2 *n 2 1. 

N n „2 Ag RI n ,f-5 ;3 ,^2 + jf2 n 2 3 _ 


= 0 


(52) 


After a particular value of kg is found from equation (52), the 
shape of the corresponding antisymmetrical free-free mode may be obtained 
by giving C the arbitrary value of unity and calculating the b's 
and a's successively from equations (5l), (50), and (49). 


Discussion of parameters .- The parameters entering in the frequency 
equations merit discussion. The unknown natural frequency is contained 
only in the frequency coefficient kg, which is defined by the formula 


03 = k- 


El 


B ^ 


and is in common use in beam- vibration analysis. 


The param- 


eters kg and k^j are identical with the shear and inertia parameters 

defined in reference 4, which considers the effect of only transverse 
shear and rotary inertia on beam vibrations. The quantity Ag which 
appears in the present definition of kg is actually the effective 
shear- carrying area when plane sections are constrained to remain plane; 
that is, when shear lag is neglected. The remaining parameters appearing 
in the present derivation, .namely, . A/Ag, K, and A]_, A . * . are 


essentially shape parameters which actually depend only on the contour 
of the cross section; as shown in the appendix. 


r a 2 

A .2 z — • n n 


*' -•.A ...2 


n=l,2,3 
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and 



and the A n ’s are simply the Fourier coefficients of the function sin 6, 
which is dependent only on the shape of the cross section. These shape 
parameters are related to shear- lag effects and their interaction with 
transverse shear and longitudinal inertia. 

The effect of longitudinal inertia is jassociated with the parameter 
k R j. If the effect of longitudinal inertia is to he neglected, it is 

sufficient to set k R j equal to zero in the final frequency equation. 

If is equal to zero, becomes independent of kg. Appreciable 

simplification in a trial-and-error solution for the natural frequency 
then results since, with Bi independent of kg, the infinite s umma tion 
contained in % is also independent of kg and need be calculated 
only once for any particular beam. As is shown in the following section, 
the effect of disregarding the influence of longitudinal inertia may 
often be negligible. 

Without presentation of details, it may be mentioned that for the 
case of a circular cylinder, which has no shear lag, all the A n ’s 
except Ap vanish and the frequency equations (to), (4l), and (52) may 
be put into closed forms identical to those given in. reference 4. Again, • 
if in the general frequency equations kg is set equal to zero, the 
equations may be put into closed forms equivalent to those of refer- 
ence 4 where only rotary inertia is considered. 


. RESULTS FOR CYLINDRICAL BEAMS OF RECTANGULAR CROSS SECTION 


In order to show quantitatively the effects of shear lag, transverse 
shear deformation, and longitudinal inertia on the natural frequencies 
of hollow thin-walled cylindrical beams, numerical calculations have been 
performed for cylinders of rectangular cross section oscillating as free- 
free beams. The calculations have been limited to symmetrical modes of 
vibration, and consequently the frequency equation (4l) is applicable. 

For rectangular cross sections the quantity N^ may be put into closed 
form as shown in the appendix, and this closed- form version of Nj. was 

used in the calculations. A' value of E/G equal to 2.65 (appropriate 
for aluminum alloys) was assumed. 
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The results of these calculations are shown in figures 2, 3 , and 4. 

In figure 2 , the ratio of the natural frequency cd to the natural fre- 
quency Ob obtained from elementary beam theory is shown as a function 
of the plan-form aspect ratio L/b for cross-sectional aspect ratios of 
1.0, 3.6, and 00 . The contribution -of each of the secondary effects to 
the total reduction in the natural frequency for the cross-sectional • 

aspect ratios - = 3.6 and 1.0 can be seen in figures 3 and 4, respec- 
a 

tively. The cross-sectional aspect ratio of 2. = co corresponds to the 

limiting case of a beam where the effects of transverse shear deformation 
and longitudinal' inertia are negligible' and therefore the reduction in 
natural frequency is due entirely to shear lag. 

The dashed lines in figures 3 and 4 show the reduction in frequency 
due to the inclusion of the effect of only transverse shear deformation 
as obtained from reference 4. 

The long- and short-dash lines are calculated from the frequency 
equation (4l) with k R j =0 and consequently represent the reduction 

in natural frequency when both shear lag and transverse shear deformation 
are taken into account. Thus the hatched area between the dashed and 
the long- and short-dash lines may be considered as showing the additional 
reduction in natural frequency when the influence of shear lag is con- 
sidered. Finally, the solid lines are calculated with k R j taken into 

account, and consequently the shaded area shows the additional influence 
of longitudinal inertia in reducing the frequency. 

Examination of figures 3 and 4 and the curves for ^ = 00 in fig- 
ure 2 shows that the influence of shear lag increases as the cross- 
sectional aspect ratio increases; whereas the influence of transverse 
shear and longitudinal inertia decreases with increasing cross-sectional 

aspect ratio. Indeed, it appears from the results for — = 3.6 that 

a 

for this aspect ratio the effects of longitudinal inertia may already 
he considered practically negligible. 

A word of caution concerning the interpretation of figures 3 and 4 
may be in order. Since in some cases the depth of the hatching increases 
with increasing L/b, it might appear, at first glance, that the shear- 
lag effect increases with increasing plan-form aspect ratio. However, 
if the additional effects of shear lag are considered on a percentage 
basis with the dashed line as a base, it will be found that shear-lag 
effects actually reduce in percentage with increasing L/b. A similar 
criterion should be used in judging the influence of longitudinal 
inertia. 



20 


NACA TN 2682 


CONCLUDING REMARKS 


The numerical calculations show that secondary effects have appre- 
ciable influence on the natural frequencies of rectangular box beams of 
uniform wall thickness. These results constitute an indication of the 
probable inadequacy of elementary beam theory for the vibration analysis 
of actual aircraft structures of the monocoque and semimonocoque type 
and emphasize the need for practical calculation procedures for such 
structures that would take into account transverse shear deformation, 
shear lag, and, when necessary, longitudinal inertia. The general solu- 
tions presented for cylinders of uniform thickness, as well as the 
numerical results for rectangular box beams, should be useful in the 
assessment of the accuracy of any procedure of this kind that may be 
developed. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. ,. January 21, 1952 
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APPENDIX 


TRANSFORMATION OF PARAMETERS 


Expressions for Aq/a, I, and K 2 .- If sin 6 is expanded into a 
Fourier series 


sin 0 


00 



n=l,2,3 


sin 


2nns 

P 


(Al) 


the Fourier coefficients A n are the same as those defined in equa- 
tion (24); that is, 


An 


2 

P 



sin 6 sin 


2n«s 

P 


ds 


(A2) 


The effective shear area Ag (equation (25)) can now he written as 
a function of the Fourier series expansion for sin 6 as 



(A3) 


With the use of the appropriate orthogonality conditions, equa- 
tion (A3) becomes, after the integration is performed. 


A S = 51 A n 2 = £ 51 A n 2 


n=l,2,3 


n=l,2,3 


or 


As 

A 


1 

2 


00 



n =l,2,3 



/ 


(A4) 



22 


NACA TN 2682 


The moment of inertia I of a cylinder is defined as (see fig. l) 



ds - ' Ay 2 


(A5) 


where y is the y-distance to the center of gravity of the cross section 
and is given by ' . 


But 


J 0 

7 = — 


-i: 


yt ds 


y = / sin 6 ds 

'o 


(A6) 


(A7) 


or 


• X,3 ’■=('- “ “) 


' (A8) 


and, consequently, 


y = 


An 


n=l,2,3 


P 

2nn 


(A9) 


With the use of equations (a 8) and (A9) , the expression for I in 
equation (A5) becomes ■ , 


I = 


kg c 


— ' A 2 

A n 


8 " 2 n=^273 n2 


(A10) 


With the series expansion for I .in equation (A10), the parameter K 2 , 
as defined in equation (29), becomes 

\ " i t 

‘ ! . 00 O 

tc 2 = -L A. A n 

jf2 Ac to r) n2 

11 n=l,2,3 


(All) 
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Transformation of expression for N^.- In equation (37) was 

defined as 


N ± = 


iV 


• 2 2 
ir*r 


00 

Z 


2 a 2 


n A 


n 


8kg 2 16 kg 2 n=lj2,3 ks^j 2 + ^n 2 


I k B 2 


(A12) 


The infinite series that appears in this expression converges as 
A n 2 and therefore is a relatively slowly converging series. In order 
to increase its rate of convergence, the following transformations are 
made. ; 

By adding and subtracting A^/k 2 inside the infinite summation in 
equation (A12) and using equation (a 4), the equation simplifies to 


Ni = 


16 Ag 


00 

rr*i 2 E 


A 2 
^n 

n=V2,3 ks 2 ®! 2 + 


J^n 2 


Ik 2 

" 2 k b 


(A13) 


By adding and subtracting A n 2 ^K 2 n 2 inside the infinite summa- 
tion in equation (A13) and using equation (All), the expression for 
can be transformed to 


N,- = 


i 2 * 4 

32 


B,- 


i 2 n 2 A 


16 A S 


4 - bA 3 


I 


A n 2 


n=lj2,3 K 2 n 2 (kB 2 Bi 2 + K 2 n 2 ) 


-iv 


( A14) 


The infinite series in equation (Al4) converges as A n 2 /n^ and 

therefore is considerably more quickly convergent than the series in equa- 
tions (A12) and (A13) , which converge as A n 2 and A n 2 /n 2 , respectively. 


Closed form of for cylindrical beams of rectangular cross 

section .- For a cylindrical beam of rectangular cross section, with 
dimensions as shown in figure 2, it is possible to write the expression 
for in a closed form. The parameters for such a cross section 

become 
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Ag = 4at 

A=4(a+b)t=pt 
An = 0 


4 2nnb 

— cos 

rut p 


(n even) f 
( n odd ) 


(A15) 


With equations (A15) the parameter shown in equation (A12) 

becomes 


Ni = 


i 2 * 2 


1 £ 


cos 


2 2nnb 


8*s 2 '“‘B 2 8 „*73,5 g? Bl 2 + n 2 


- 5 k B 2 


(A16) 


or 


N ± = 


i 2 * 2 i 2 


/ 


r- 00 — 


51 


8kg 2 8kg 2 S \ n= ^ Ts kg 2 


\n=l,3,5 -S_ b ± 2 + n 2 
' K 2 


00 

z 


cos 


4mtb 


n=l>3,5 =§- Bi 2 + n 2 
K 2 


. 1 k 2 
2 B 


(AIT) 


Each of the infinite summations in equation (A17) can now be written 
in closed form as shown in reference 5 , and the closed expression for % 
then becomes 


N, = 


i 2 rt 

8kg 2 


It - 


K p 
4kgB i a 


sinh j j. B ± (& - l) 
c °sh f ^ B i 


tanh 




(A18) 
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(c) Cantilever beam. 


(d) Free-free beam. 


Figure I. Coordinate systems and sign conventions 




Figure 2— Change in the natural frequency of a symmetrically vibrating free-free cylinder due to the inclusion of secondary 

effects. 










